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Abstract 


This report presents a theoretical analysis of the dynamic effects of angular 
motion in yaw and pitch on model attitude measurements in which inertial sensors 
were used during wind tunnel tests. A technique is developed to reduce the error 
caused by these effects. The analysis shows that a 20-to-l reduction in model attitude 
measurement error caused by angular motion is possible with this technique. 


Introduction 

The standard model attitude measurement instru- 
mentation used at Langley Research Center (LaRC) is 
the inertial sensor or angle-of-attack (AOA) package 
described in detail in reference 1. Installed in a model 
during wind tunnel tests, the AOA package generates an 
output signal proportional to the model pitch attitude. 
This signal is then heavily low-pass filtered before final 
data processing to remove time-varying components 
caused by random flow noise and vibration. Certain test 
conditions produce oscillatory model and sting motion in 
pitch and yaw, which produces a centrifugal or normal 
acceleration component in which the average value 
results in constant offset errors in the sensed AOA mea- 
surement. The magnitude of the offset can be as great as 
0.1°, which is unacceptable in AOA measurements in 
which the accuracy requirement can be as low' as ±0.01 
This paper presents a dynamic analysis and proposes a 
spectral analysis technique for detecting model and sting 
motion and for correcting the resulting AOA errors by a 
factor of 20 to 1 . 

Symbols 

Ap pitch oscillation amplitude, rad 

A yaw' oscillation amplitude, rad 

a normal acceleration component caused by 

pitch motion, in/sec~ 

a normal acceleration component caused by 

yaw motion, in/sec^ 

E p pitch plane accelerometer 

E standard AOA sensor 

r 

E yaw plane accelerometer 

g gravitational constant, in/sec 2 

J ( u) nth-order Bessel function of first kind 
r pitch motion radius, in. 

r yaw motion radius, in. 

S AOA sensor sensitivity, V/g 

t time, sec 

u parametric constant, — 


V aoa corrected, filtered sensor voltage output, V 

Vppj power spectral line amplitude, 

frequency 20 , V 

V fl | uncorrected, filtered sensor voltage output, V 

l unf uncorrected, unfiltered sensor voltage 
output, V 

angular pitch plane position, rad 

P angular yaw' plane position, rad 

9 model pitch attitude, rad or deg 

a A ~ A standard deviation of corrected, filtered sensor 
AOA 

output, deg 

Oppj. standard deviation of power spectral line, 
frequency 20 , deg 

Gpp-p . standard deviation of power spectral line, / th 
oscillatory mode, deg 

O f j standard deviation of uncorrected, filtered 
sensor output, deg 

longitudinal oscillation frequency, rad/sec 
0^ pitch plane oscillation frequency, rad/sec 

0 yaw plane oscillation frequency, rad/sec 

0) angular velocity in pitch, rad/sec 

0) angular velocity in yaw, rad/sec 

Abbreviations: 

AOA angle of attack 

FFT fast-Fourier transform 

Analysis of Yaw and Pitch Motion 

As shown in figure 1, let \$ p and P v denote the 
angular motion of the model in pitch and yaw, respec- 
tively. Let 0 p and 0 V , and A p and A y denote the fre- 
quencies and amplitudes of oscillation in the respective 
planes. Yaw and pitch will be analyzed separately. 



Motion in Yaw Plane 

Yaw motion is first modeled as in the following 
equations: 

B = A sin6 t i 

.v v Y v 

co y = 4 v $ y cos0 v r l (1) 


where co y is the angular velocity, a ny is the normal com- 
ponent of the acceleration vector, r' y is the radius, and t 
is time. Combining the second and third expressions of 
equations (1) yields the normal component of the accel- 
eration vector as a function of time: 
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a ny (0 = V 2 V ( 1 + cos2( l y) (2) 

The standard LaRC AO A package is installed in an 
airplane model as shown in figure 2, with its sensitive 
axis parallel to the longitudinal axis of the model. The 
AOA output signal caused by angular motion in yaw, 
prior to filtering, is 
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g sin0 
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( 1 + cos2 4> v r) 


(3) 


where V unf is the uncorrected, unfiltered sensor-package 
output voltage, S is the AOA sensor sensitivity, g is the 
gravitational constant, and 0 is the model pitch attitude 
or angle of attack, expressed in radians for computation 
and converted to degrees for final presentation. The three 
terms of equation (3) are presented in figure 3 as func- 
tions of time. After filtering, the sensor output is the sum 
of the first two terms of equation (3) as shown in the fol- 
lowing equation: 


v a = s 
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where is the uncorrected, filtered sensor output. The 
offset error, given by the second term r y A y ^ y /2 and 
equal to the amplitude of the third term of equation (3), 
can be determined from the power spectrum of the unfil- 
tered sensor output at frequency 2<j) v . Thus, 


V = S 

V ppp o 


( A 2.2 
r A*v 


(5) 


where Vppj is the power spectral amplitude at frequency 
2<J) V . The corrected, filtered sensor output V AOA finally 
becomes 


V = V + V 

V AOA V fi\^ V FFT 


(6) 


Motion in Pitch Plane 

Pitch motion is similarly modeled in the following 
equations: 


e, 
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= A sincb t 
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= A 6 cos<b t l 
p r y y p 
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where to^ is the angular velocity, a np is the normal 
acceleration component, r p is the radius, and t is time. 
The normal acceleration component can be expressed as 


up 


2 2 
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If the angular pitch motion is superimposed on the 
pitch angle, equation (3) becomes 
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Equation (9) can be further approximated using 
equation (6) to yield the following: 
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In addition to the normal acceleration vector component, 
comparison of equations (3) and (10) shows the presence 
of another sinusoidal component at frequency caused 
by frequency modulation. The modulation component at 
<)>p does not produce offset, and error correction is 
unnecessary at this frequency. Additional sensors may be 
required to identify these modulation frequencies. The 
same error correction method is used for both yaw and 
pitch. That is, the corrected AOA sensor output is the 
sum of the uncorrected filtered output and the amplitude 
of the spectral component at frequency 2§ p . 


Magnitude Estimate of Error Caused by Angular 
Motion 

The AOA offset error, as expressed in equation (5), 
is presented in figure 4 as families of curves parameter- 
ized for typical values of the angular motion amplitude 
and its radius over frequencies from 0 to 30 Hz. Experi- 
mental AOA offset errors of 0.1° have been observed. 
These errors correspond to a 0.2° motion amplitude at a 
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30-in. radius at 10 Hz. Figure 4 illustrates curves of AO A 
offset errors at radii of 6 and 30 in. and at angular ampli- 
tudes from 0.1° to 0.2°. Other combinations of radius, 
amplitude, and frequency may produce AOA offset 
errors up to 0.1 °. 

AOA Sensor-Package Configuration 

The total sensor output which results from multiple 
oscillatory modes is presented and shows that additional 
sensors are required to identify uniquely those oscillatory 
modes arising from angular motion. An error analysis 
follows that estimates the error reduction provided by the 
proposed technique. 

AOA Sensor- Package Output 

The AOA sensor-package output contains at least the 
following signal components: 

1. Constant component sin0 caused by model pitch 
angle 

2. Time-varying and constant components caused by 
angular yaw motion 

3. Time-varying and constant components caused by 
angular pitch motion 

4. Time-varying components along the model longitu- 
dinal axis 

Note that all time-varying components can be removed 
by low-pass filtering. Spectral analysis of the unfiltered 
sensor output identifies all time-varying components in 
the signal. Appropriate correction is applied to those 
components producing AOA error. Term-by-term error 
correction is required for multiple oscillatory modes in 
pitch and yaw identified by the spectral analysis. 

Three- Axis Sensor Package 

The proposed technique to reduce AOA error 
requires a spectral line estimate of the unfiltered AOA 
sensor-package output. However, as indicated in 
equation (10), indiscriminate use of fast-Fourier trans- 
form (FFT) spectral lines may result in improper correc- 
tion for modulation components. Consequently, 
additional sensors in pitch and yaw are required to iden- 
tify only those oscillatory components producing AOA 
errors. Figure 5 shows a proposed three-axis sensor pack- 
age that consists of a standard LaRC AOA sensor E r 
aligned with the model longitudinal axis and two minia- 
ture accelerometers £ v and E p mounted in the yaw' and 
pitch planes, respectively. Accelerometer £ v , mounted 
transversely to the roll axis, responds only to yaw 
motion. Since £ v does not respond to normal accelera- 


tion produced by angular motion in either pitch or yaw, 
its output does not contain the constant terms or the time- 
varying components at frequencies 2<f> v and 2§ p appear- 
ing in equations (2) and (8). Accelerometer £ v output, 
however, does contain the tangential acceleration com- 
ponent at frequency 0 V . Similarly, accelerometer E p , 
mounted orthogonally to £ r , does not respond to normal 
acceleration produced by either pitch or yaw motion. It 
does respond to the tangential acceleration component in 
the pitch plane at frequency § p and to the cosine compo- 
nent of the gravitational vector. The presence of angular 
motion in each plane can thus be identified from the 
power spectral component frequency signatures of the 
three accelerometer outputs, as summarized in table I. 


Table I. Sensor Power Spectral Frequency Components 


Input motion 

Components for accelerometer — j 

E r 

fv 

E 

p 

Axial 

«•/ 



Yaw 

2*v 

4> v 


Pitch 

V 2d V 


0 


Of the four spectral components in the £,- power spec- 
trum, only those modes at frequencies 2<J> V and 2§ p pro- 
duce AOA offset errors that require correction. 
Therefore, the noncontributing modes at <(>/ and § p must 
be identified and discarded. The following procedure 
identifies those frequency components in E r to be 
included or eliminated for correction of AOA offset 
errors. 

1. Frequency ty p is obtained from the cross spectrum 
of E p and E r outputs; <b p is eliminated, and a cor- 
rection is made at 2§ p . 

2. Frequency <j> v is obtained from the £ v power spec- 
trum. The E r power spectrum is searched for a 
component at frequency 2<j> v for correction. 

3. The remaining frequency components in the E r 
spectrum, assumed to be longitudinal, are 
eliminated. 

Measurement Uncertainty Estimate 

The variance of AOA measurements provided by 
this method is obtained from the variance of 
equation (6): 

a AOA “ a fil + CT FFT (11) 
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If additional oscillatory modes that contribute signifi- 
cantly to the total AOA error are identified, their individ- 
ual error values must be included in the following form: 

n 

a FFT = X a FFT,/ U2) 

/ = 1 

2 

where Oppp is the measurement error variance of the 
zth oscillatory mode. 

For a maximum AOA offset error of 0.1° and an 
estimated 5-percent spectral analysis accuracy, the stan- 
dard deviation of the AOA offset error becomes 

= 0.05x0.1° = 0.005° (13) 

Thus, a 20-to-l reduction in the error portion caused 
by angular motion is possible. Equation (11) and the 
assumption that = ±0.01° (ref. 1) show that the 
total error is 

"AOA = °- OII ° " 4 > 


Concluding Remarks 

This paper is a theoretical treatment of angular 
motion effects upon inertially sensed model attitude mea- 
surements. A technique was developed to identify the 
time-varying components that produce angle-of-attack 
(AOA) offset errors. The analysis quantifies the offset 
error magnitudes, based on knowledge of the angular 
motion radii and frequencies. Spectral analysis tech- 
niques identify only those frequency components of the 
motion that contribute to AOA offset errors. A technique 
that does not require prior knowledge of the radius of 
motion is proposed to estimate AOA offset error magni- 
tudes. Experimental verification of the proposed method 
has not been attempted. Note that the AOA offset error 
cannot be removed by filtering or time averaging. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
November 28, 1994 
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Appendix 

Approximation of Angular Motion in Pitch 

Equation (10) describes the angular pitch motion. 
The first term contains sin0 , and is phase modulated by 
the sinusoidally varying angular motion at amplitude 
Ap/r p at frequency § p , as shown in the following 
equation: 


f(t) = sin ( 0 + u sin<j |> /) (Al) 

where u = A p /r p . Expand equation (Al) using trigono- 
metric identities into the form 

f(t) = sin0cos (ttsin<() /) + cosOsin (wsincjw) (A2) 

The second factor of each term of equation (A2) can be 
expanded as follows: 


cos (i/simjw) = J 0 (u) + 2 ^ J 2n (u) cos2 § p t 

n = i 

sin (Hsin<|> r) = 2 ^ J 2n , (a) sin (2n - 1 ) § p t 

n = 1 


(A3) 


where J n (w) is the /?th-order Bessel function of the first 
kind. For the values of A p and r p , cited in figure 4, 
u ranges from 0.001 to 0.008, for which the Bessel func- 
tions are closely approximated as 


y 0 (w) * l , 

y, («> * | 

J„(U) =0 (»>D J 


(A4) 


where 7 0 (0.001) = 1.00000. 7 0 (0.008) = 0.99998, 
7,(0.001) = 0.00050000 , 7,(0.008) = 0.0040000. 

J 2 (0.001) = 1.2x10 ' , J 2 (0.008) = 0.0000008, 

and so on. Therefore, equations (A3) can be approxi- 
mated as 


cos (u sin 0/) «J n (w) ~ 1 1 

(A5) 

sin (w sinter) « 2J j (u) sin 0^ = wsinb/ J 

Thus, function f(t) in equation (A2) simplifies to 

f{t) ~ sin 0 + w cos Osin (A6) 
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Figure 3. AOA sensor output. 







REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden tor this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathenng and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect ot this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services. Directorate for Information Operations and Reports. 1215 Jefferson 
Davis Highway, Suite 1204. Arlington, VA 22202-4302, and to the Office of Management and Budget. Paperwork Reduction Project (0704-0186), Washington, DC 20503. 

1 AGENCY USE ONLY (Leave blank) |2. REPORT DATE 1 3. REPORT TYPE AND DATES COVERED 


February 1995 Technical Memorandum 


4. TITLE AND SUBTITLE 

Effects of Yaw and Pitch Motion on Model Attitude Measurements 

5. FUNDING NUMBERS 

WU 505-59-54-02 

6. AUTHOR(S) 

Ping Tcheng, John S. Tripp, and Tom D. Finley 


7. PERFORMING ORGANZATION NAME(S) AND ADDRESS(ES) ~~ 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 

NASA Langlev Research Center 
Hampton, VA 23681-0001 

L-17381 

9. SPONSORIING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 

NASA TM-4641 

11. SUPPLEMENTARY NOTES 

12a. DISTRIBUTION/ AVAILABILITY STATEMENT 

Unclassified-Unlimited 
Subject Category 35 

Availability: NASA CAST (301) 621-0390 

12b. DISTRIBUTION CODE 


13. ABSTRACT (Maximum 200 words) 


This report presents a theoretical analysis of the dynamic effects of angular motion in yaw and pitch on model atti- 
tude measurements in which inertial sensors were used during wind tunnel tests. A technique is developed to 
reduce the error caused by these effects. The analysis shows that a 20-to-l reduction in model attitude measurement 
error caused by angular motion is possible with this technique. 


14. SUBJECT TERMS 15. NUMBER OF PAGES 

Model attitude measurement; Inertial angle-of-attack sensors; Dynamic attitude 9 

measurements 16. price code 

A02 

17. SECURITY CLASSIFICATION 1 8. SECURITY CLASSIFICATION 1 9. SECURITY CLASSIFICATION 20. LIMITATION 
OF REPORT OF THIS PAGE OF ABSTRACT OF ABSTRACT 

Unclassified Unclassified Unclassified 


NSN 7540-01-280-5500 


Standard Form 298 (Rev. 2-89) 

Prescribed by ANSI Std. Z39-18 
298-102 

























